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Abstract. The primordial abundance of deuterium produced during Big Bang Nucleosynthesis 
(BBN) depends sensitively on the universal ratio of baryons to photons, an important cosmological 
parameter probed independently by the Cosmic Microwave Background (CMB) radiation. Obser- 
vations of deuterium in high-redshift, low-metallicity QSO Absorption Line Systems (QSOALS) 
provide a key baryometer, determining the baryon abundance at the time of BBN to a precision of 
~ 5%. Alternatively, if the CMB-determined baryon to photon ratio is used in the BBN calcula- 
tion of the primordial abundances, the BBN-predicted deuterium abundance may be compared with 
the primordial value inferred from the QSOALS, testing the standard cosmological model. In the 
post-BBN universe, as gas is cycled through stars, deuterium is only destroyed so that its abundance 
measured anytime, anywhere in the Universe, bounds the primordial abundance from below. Con- 
straints on models of post-BBN Galactic chemical evolution follow from a comparison of the relic 
deuterium abundance with the FUSE-inferred deuterium abundances in the chemically enriched, 
stellar processed material of the local ISM. 
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INTRODUCTION 

Of the light nuclides synthesized in astrophysically interesting abundances (D, 3 He, 
He, 7 Li) during Big Bang Nucleosynthesis (BBN), the post-BBN evolution of deu- 
terium is the simplest. As gas is cycled through stars, deuterium is only destroyed [Q]], 
so that its abundance, observed anywhere in the Universe, at any time in its evolu- 
tion, is constrained to be bounded from above by its primordial, BBN abundance: 
(D/H)obs < (D/H)p. In systems at high redshift and/or with very low metallicity, the 
observed deuterium abundance should approach its primordial value. Its simple post- 
BBN evolution, along with the sensitivity of its predicted BBN abundance to the baryon 
to photon ratio % = n^/n r ((D/H)p oc r/ B L6 ), identifies deuterium as the baryometer of 
choice. The baryon mass density parameter, ^b (the fraction of the present critical mass 
density contributed by baryons), is related to t]b and the present value of the Hubble 
parameter, H = 100/z kms^Mpc -1 , by Q.^,h 2 = T] 10 /274, where T]io = 10 10 7] B - Given 
the dependence of (D/H)p on 7]b, observations which constrain the relic deuterium abun- 
dance to, say, ~ 10%, lead to a ~ 6% measurement of the baryon to photon ratio when 
the Universe was only a few minutes old. Alternatively, the baryon to photon ratio de- 
termined by observations of the cosmic microwave background (CMB), which probe a 



FIGURE 1. Deuterium abundances, the D to H ratios by number, derived from observations of high 
redshift, low metallicity QSO Absorption Line Systems, are shown as a function of the corresponding 
neutral hydrogen, H i, column densities. The solid line shows the weighted mean of the D/H ratios and the 
dashed lines indicate the ±1(7 error estimates; see the text. 



time when the Universe was some 400 thousand years old, may be used in the BBN 
calculation to predict the primordial D abundance, which may then be compared to the 
relic value inferred from observations of deuterium at high redshifts in systems of very 
low metallicity, testing the standard model of cosmology. 

In the post-BBN Universe, the D "astration" factor, /d = (D/H)p/(D/H), measures 
the virgin fraction (1//d); 1 / Id is the fraction of gas which has never been processed 
through stars. Due to stellar nucleosynthesis, as the metallicity, Z, in a system increases, 
the deuterium abundance should decrease, suggesting that, in the absence of dust deple- 
tion, the observed values of D/H and Z should be anti-correlated. 

BBN AND THE PRIMORDIAL ABUNDANCE OF DEUTERIUM 

In Figure [TJ are shown the deuterium abundances inferred from observations of seven, 
high redshift, low metallicity QSO Absorption Line Systems (QSOALS) [2]. The 
weighted mean of these abundance determinations provides an estimate of the primor- 
dial deuterium abundance, vdp = 10 5 (D/H)p = 2.7 ± 0.2. Since the dispersion of the 
individual abundance determinations around the mean value is very large, the formal 
error in the mean has been multiplied by the square root of the reduced x 2 to provide a 
more realistic error estimate. 

From BBN |Ht], for yrjp = 2.7 ± 0.2, the baryon abundance is determined, at the ~ 5% 
level, to be 7710 = 6.0 ±0.3. This determination, some ~ 20 minutes after the expansion 
has begun, is in excellent agreement with the value determined independently from the 
CMB, some 400 kyr later, rjio = 6.1 ±0.2 (see Simha & Steigman [3] and references 
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FIGURE 2. The FUSE-observed ISM abundances of deuterium versus those of oxygen observed along 
the same lines of sight; from 0]. 



therein). If the CMB-determined value of the baryon abundance is used to predict the 
primordial abundance of deuterium, the result is vdp = 2.6 ±0.1, in perfect agreement, 
within the estimated errors, with the relic abundance inferred from the QSOALS data. 

OBSERVATIONS OF D IN THE ISM 

The conventional wisdom (ante?) has been that mixing of the ISM within 1-2 kpc 
of the location of the solar system is efficient so that this relatively local gas is well- 
mixed with homogenized abundances. In particular, the conventional wisdom led to the 
expectation that observations of deuterium and, e.g., oxygen within the local (~ 1 kpc) 
ISM would reveal uniform abundances, "the interstellar abundances", of these (and 
other) elements. Beginning with the pioneering Copernicus observations [|4J], continuing 
with observations provided by the IUE satellite, and the GHRS and STIS instruments 
onboard the HST, and confirmed by those of the FUSE mission |50, the data reveal that 
nothing could be further from the truth. As may be seen in Figure [2l the deuterium 
abundances inferred from observations in the relatively local ISM (see [50 for details 
and further references) vary over a range of a factor of ~ 3, while those of oxygen range 
over a factor of ~ 4. 

While there has not been much discussion in the literature of the factor ~ 4 range in 
oxygen abundances observed along lines of sight (LOS) within ~ 1 — 2 kpc of the Sun, 
much attention has been paid to the factor of ~ 3 range in the deuterium abundances. 
The general consensus is that the large range in observationally-inferred deuterium 
abundances is the result of preferential depletion onto grains of D, relative to H JfJ. 
Although preferential depletion of D is likely the correct explanation for most, possibly 
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FIGURE 3. In the left-hand panel the FUSE observed D/H ratios are shown versus the corresponding 
LOS values of E(B - V) and, in the right-hand panel of log N(H 2 ); from 0]. 



all, of the observed variation in ISM D abundances, there is, in fact, surprisingly little 
independent data in support of it. For example, if depletion onto grains is the culprit, 
it might be expected that the D/H ratios would correlate (or anticorrelate) with the 
reddening or, with the amount of molecular hydrogen. However, as the two panels of 
Figure [3] reveal, both high and low D abundances are found at very small values of 
E(B - V) and of N(H2) and, while there may be weak evidence in favor of decreasing 
D abundances at high E(B - V), both high and low D/H values are found along the 
LOS with higher N(H2). This absence of a smoking gun for depletion is also revealed 
in Figure @] where the D abundances are shown as a function of the corresponding iron 
abundances along the same LOS. While there IS evidence of a correlation between D/H 
and Fe/H at low Fe abundances, that trend disappears for log ype = 10 6 (Fe/H) > 0, 
hinting that along those LOS with log ype ~ 0, deuterium may be undepleted. 

Whatever is the correct explanation for the large range in the ISM abundances of D 
and O, it is clear that the local ISM is NOT homogeneous. Since any local depletion of D 
onto grains has not been homogenized in the ISM, it is not unreasonable to suspect that 
there might also be some LOS along which infalling, D-enhanced, O-poor gas has been 
incompletely mixed with the processed, D-astrated, O-enhanced gas already present 0]. 
This latter effect is illustrated in Figure 0] by the solid curve which shows the effect of 
unmixed, nearly primordial infall [7] to yield higher D and lower Fe abundances than 
is predicted by the fiducial Galactic chemical evolution (GCE) model, shown by the 
star. Abundances to the left and below the dashed lines may have resulted from Fe and 
D depletion respectively, while any D abundances above the dashed line may have been 
contaminated by infall. These two competing processes, depletion onto grains and infall, 
complicate the attempt to use the data to infer the "true" ISM abundances of D and Fe 
(and O). 
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FIGURE 4. FUSE observed ISM D and Fe abundances (y D = 10 5 (D/H), y Fe = 10 6 (Fe/H)); from @]. 
The star shows the ISM abundances predicted by the fiducial GCE model (see the text). D and Fe 
abundances below and to the left of the dashed lines are consistent with depletion. The solid curve shows 
the effect of local abundance fluctuations resulting from incompletely mixed infall of D-enhanced, Fe-free 
material. 



GALACTIC CHEMICAL EVOLUTION 

Unable to survive the high temperatures in stars, deuterium is destroyed when gas is 
incorporated in stars. As gas cycles through stars and the stellar processed material 
is returned to the interstellar medium (ISM) of the Galaxy, the metallicity of the ISM 
increases while the abundance of deuterium decreases. Infall to the disk of the galaxy 
of essentially unprocessed gas plays a crucial role in the chemical evolution of the 
Galaxy [|80. Infall of primordial (or very nearly primordial) gas which is metal-free (or 
very nearly so) and whose deuterium abundance is primordial (or very nearly so) dilutes 
the metallicity of the ISM while enhancing its D abundance 10, & Gil- 

In Steigman, Romano, Tosi (SRT) ||7j], we considered GCE models with several dif- 
ferent prescriptions for the stellar IMF and stellar lifetimes and predicted the ISM D as- 
tration factor and the current ISM oxygen (and iron) abundances. For those GCE models 
consistent with other independent, observational data, SRT found that 1.4 < /d 1.8 
and that 7 <y ^ 12, where yo = 10 4 (O/H). For our fiducial model, SRT adopted the 
Scalo IMF ffllTl along with the Schaller etal. lfl2ll prescription for stellar lifetimes, corre- 
sponding to /d = 1-39 and yo = 8.4; these values are indicated by the stars in Figures |4] 
and|5] As noted by SRT if recent infall of unprocessed material were not fully mixed in 
the ISM along some LOS, there would be an anticorrelation between the D and O abun- 
dances along such contaminated LOS, predicting Vp BS = 2.61 — 0.087vq BS (the dashed 
line in the upper panel of Fig. [5]). SRT also noted that since some of the dispersion ob- 
served among the D and O abundances may be due to systematic errors in the H I column 
densities, the D/O ratio would be unaffected by such errors. For LOS contaminated by 




FIGURE 5. The upper panel shows the FUSE observed deuterium and oxygen abundances, while the 
lower panel plots the D to O ratios as functions of the oxygen abundances; from |7|]. The stars show the 
ISM D and D/O versus O abundances predicted by the fiducial GCE model (see the text). Abundances 
below and to the left of the dashed lines are consistent with D and O depletion. 



unmixed infall, 10(D/O) = vd/vo, 10(D/O) = 26. 1 - 0.87 (the dashed line in the lower 
panel of Fig. O). 

TOWARDS AN ESTIMATE OF THE D ASTRATION FACTOR 

Estimating the Galactic D astration factor is complicated by the possibilities that the 
"true" ISM D abundance along some LOS may have been reduced by preferential D 
depletion onto grains and/or, enhanced by incompletely mixed infall of gas with the 
higher, primordial D abundance. Linsky et al. Q5D assume that unmixed infall makes 
no contribution at all so that, allowing for depletion, the maximum of the ISM D/H 
ratios provides a lower limit to the "true" ISM D abundance. Linsky et al. then find 
the weighted mean of D/H for the 5 highest D abundances, concluding that y^ M > 
yMAX _ 2. 17 ±0.017 (or, including a correction for the Local Bubble D abundance, 
-yMAX _ 2 37 ±0.024). These choices correspond to astration factors /d < 1.1 — 1.2, 
considerably smaller than those predicted by the GCE models, jfo « 1.4—1.8 |0]. 
In contrast, SRT allow for the possibility that the highest observed D/H ratios may 
represent the "true" ISM D abundance, although unmixed infall may have contaminated 
the D abundance along some LOS and, they note that there are many more than 5 LOS 
with high D/H abundances which are equal within the errors. SRT find that for the 18 
FUSE LOS with the highest D/H ratios (vrj > 1.5), nine have vd > 1.9 and nine have 
< 1 - 7, and they find for the weighted mean, M = 1.9. For the primordial abundance 
yop = 2.7 ±0.2, this corresponds to fa = 1.4 ±0.1, in agreement with the fiducial model. 



FIGURE 6. The 68%, 95% and 99% contours in the y D - w plane (y D = 10 5 (D/H) ISM ) for a top-hat 
prior distribution for three subsets of the FUSE D/H data; see the text for details. From [13]. 
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FIGURE 7. The likelihood distribution for the ISM D abundance ( vd.ism — 1 .95^) inferred from the 
data for 38 LOS (see Fig. [6j, using the top-hat prior for w. From [13]. 



Recently, B. Fields, T. Prodanovic, and I [13] have explored a different path to 



estimating y^ M , using a Bayesian approach pioneered by Hogan, Olive, and Scully ||14fl . 



In our analysis we neglect unmixed infall and assume that the scatter among the FUSE- 
observed D/H ratios is entirely due to preferential depletion of D onto grains, so that 



our estimate likely provides an upper bound to the true value of y£ . We assume 
that the gas phase D abundances (including uncertainties) are equal to the differences 
between the true ISM D abundance and corrections, w = y^ M — v§ BS , for depletion. We 
adopt a prior for the distribution of w and maximize the likelihood to find the observed 
abundances as a function of y^ M and w. In Figure [6] are shown the contours in the y^ M 
- w plane for a top-hat distribution for w for three subsets of the data. As the number of 
data points is increased, the contours shrink and w > is strongly favored, arguing in 
support of depletion. The likelihood distribution for y^ M = 1.95+QQ9, derived from the 
top-hat prior for D/H along 38 LOS II 1311 . is shown in Figure[7J Combining this estimate 
for (a lower bound to) M with the weighted mean of the QSOALS D/H ratios [2] from 

Simha & Steigman y|], yop = 2.701q 20' resu l ts in a D astration factor /d ^ 1.38^' J 
which, within the uncertainties, is consistent with all models of GCE. 
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